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1 Abstract 

Recent observations of microlensing events in the Large Magel- 
| lanic Cloud by the MACHO and EROS collaborations suggest 

that an important fraction of the galactic halo is in the form of 
Massive Halo Objects (MHO) with mass ~ O.1M . We outline 
^ . a scenario in which dark clusters of MHO with mass ~ O.1M 

and H 2 molecular clouds form in the halo at galactocentric dis- 
tances larger than ~ 10 — 20 kpc, provided baryons are a major 
constituent of the halo. Possible signatures of our picture are 
discussed. 



September 1994 

* Supported by the Swiss National Science Foundation. 

# Partially supported by Agenzia Spaziale Italiana. 



One of the most important problems in astrophysics concerns the nature 
of the dark matter in galactic halos, as suggested by the observed flat rotation 
curves. Although various exotic dark matter candidates have been proposed, 
present limits coming from primordial nucleosynthesis allow a halo made of 
ordinary baryonic matter that should be in the form of Massive Halo Objects 
(MHO) with masses in the range 10" 7 < M/M & < 10" 1 0. Paczyhski 
suggested to detect MHO using the gravitational lens effect Q|. 

Recently, the French collaboration EROS and the American- Australian 
collaboration MACHO |§] reported the possible detection of altogether six 
microlensing events, discovered by monitoring over several years millions of 
stars in the Large Magellanic Cloud (LMC), whereas the Polish- American 
collaboration OGLE || and also the MACHO collaboration have found al- 
together at least forty microlensing events by monitoring stars located in the 
galactic bulge. Taking these results at face value an average mass ~ O.1M 
U for the MHO in the halo has been derived. 

At present the few events which have been found so far by monitoring 
the LMC do not allow to make a precise estimate of the fraction of dark halo 
matter in form of MHO, nor to infer whether the MHOs are located in the 
halo or rather in the LMC itself [[jj or in a thick disk in our galaxy [§]. It 
has also been proposed |§ that microlensing by stars within the LMC itself 
can account for the observed events. 

Our aim is to present a scenario in which the halo of the Galaxy substan- 
tially consists of dark clusters of MHO and/or H 2 molecular clouds, provided 
baryons are a major constituent of the halo. We suggest that the dark clus- 
ters form at galactocentric distances larger than ~ 10 - 20 kpc. 
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The proposed picture relies on the theory for the origin of stellar globular 
clusters advocated by Fall and Rees |TIJ and on the suggestion of Palla, 
Salpeter & Stahler |Tl[ that the lower bound on the fragment masses in a 
collapsing, metal poor cloud can be as low as 1O _2 M . Here we outline the 
main facts and for more details see ref. [12| . 



As a preliminary step, let us briefly recall the key-points of the present 
wisdom about the formation of stellar globular clusters. After its initial col- 
lapse, the proto Galaxy (PG) is expected to be shock heated to its virial 
temperature T e ~ 10 6 K. This temperature lies near a very unstable re- 



gion of the cooling curve [13], so that density enhancements should rapidly 
grow as the gas cools to lower temperatures. Fall and Rees |L0|] argued that 
overdense regions in the PG cool more rapidly than average (because the 
cooling rate by hydrogen recombination is proportional to n 2 H , where n# 
is the number density of hydrogen) and then a two phase medium forms 
with cool proto globular cluster (PGC) clouds in pressure equilibrium with 
the external hot diffuse gas. The PGC clouds, which have temperature 
T c ~ 10 4 K, are gravitationally unstable when their mass exceeds the Jeans 
mass Mj = 1.18 {k B T c /fi c ) 2 G~ 3 ^P~^ 2 , where fi c ~ 1.22 m p is the mean 
molecular weight of the primordial neutral gas and P e is the pressure of 
the fully ionized plasma that fills the PG. Typical density of the PG is 
p e ~ 1.7 x 1CT 24 (R/kpc)^ 1 g cmT ? % so that the resulting PGC cloud mass 
and radius are M c ~ 10 6 (R/kpc) 1 / 2 Mq and r c ~ (19 pc) (R/kpcj^A- 1 , 
where R is the galactocentric distance and A ~ 10 is a factor which takes 
into account the shrinking of the PGC cloud due to dissipation and collapse 
during its subsequent evolution (for more details see p^| ). The main coolants 
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below 10 4 K are molecular hydrogen and any heavy element produced by a 
first generation of stars. Molecular hydrogen, however, would be dissociated 
by various sources of UV radiation such as an active galactic nucleus (AGN) 
and/or a population of massive young stars in the center of the PG. In the 
early more chaotic phase of the evolution of the PG, in fact, an AGN is 
possibly formed at its center and/or a first generation of massive stars due 
to the disruption of central PGC clouds. The last circumstance is realized 
because, in the center of the PG, the cloud collision time is shorter than the 
corresponding cooling time. In a metal poor, dust free protogalactic gas the 
molecular hydrogen abundance fn 2 — n H 2 / n H is determined by considering 
a set of atomic and molecular processes which involve the creation of inter- 
mediate H~ and H% [jrj|. With the knowledge of fn 2 one can compute the 



cooling rate A c and the cooling time 

3p c k B T c 

Tco ° l ~ 2/i c (A c - r) ' [i) 
where the heating rate T due to external radiation sources has also been 
taken into account. The subsequent evolution of the PGC cloud depends on 
the ratio between r coo i and the gravitational infall time T// = (Sn/32Gp c ) 1 ^ 2 
that, since p c ~ 400p e Q, results to be r// ~ 1.7 x 10 6 (R/kpc) l/2 years. If 
T cooi < ^~ T ff the PGC clouds rapidly cool to a temperature T c ~ 100 K before 
the gravitational instability sets in, while for t coo \ < r// the PGC clouds cool 
and contract at the same time. On the other hand, when t coo \ > Tff the 
PGC clouds do not cool, remaining at a fixed temperature. The parameters 
that mainly discriminate between the different situations are T and A c . 

In the inner halo, because of the presence of an AGN and/or a first 
population of massive stars at the center of the PG, H 2 formation and cooling 



are heavily suppressed or delayed for a wide range of both external UV fluxes 
and/or PGC cloud densities [pj]. For this case, in which r coo i > Tff, the PGC 
clouds remain at T c ~ 10 4 K for a long time. This results in an imprinting of 
a characteristic mass of M c ~ 10 6 M . Moreover, during the permanence of 
the PGC clouds for a long time in quasi hydrostatic equilibrium, propagation 
of sound waves erases all large scale perturbations leaving only those on small 
scale. After enough H 2 has formed (/# 2 ~ 1CT 3 ), the temperature suddenly 
drops well below 10 4 K because now T cooi <C Tff. The subsequent evolution of 
the PGC clouds goes on with a rapid growth of the small scale perturbations 
that leads directly (in one step) to the formation of stars within a narrow 
mass range fi~7fl . Thus, this scenario would explain the formation of stellar 
globular clusters which are observed today especially in the inner part of the 
galactic halo. The formation of the PGC clouds could have been delayed 
until after the Galaxy was enriched by metals due to a first generation of 
stars, in this way explaining the observed absence of globular clusters with 
only primordial metal abundances and the radial gradient of metallicity in 
the galactic halo. 

In the outermost part of the halo, where the incoming UV radiation flux is 
suppressed due to the distance (so that r coo i < Tff), the PGC clouds cool more 
gradually below 10 4 K. Then, cooling and collapse occur simultaneously and 
their evolution proceeds according to the scenario proposed by Palla et al. 
[PH , leading to a subsequent fragmentation into smaller clouds that remain 
optically thin until the minimum value of the Jeans mass (< O.1M ) is 
attained. In this case we expect that, because the PGC clouds do not remain 
at a fixed temperature, there should be no imprinting of a characteristic 



4 



mass. In fact, when a PGC cloud is in a quiet ambient as at the edge of 
the PG the collapse proceeds with a monotonic decrease of the Jeans mass 
and a subsequent fragmentation into clouds with lower and lower masses. 
This process stops when the fragments become optically thick to their own 
line emission. This happens because in a metal poor cloud virtually all the 
hydrogen gas is converted into molecular form due to three-body reactions 
(H + H + H -> H 2 + H and H + H + H 2 -> 2H 2 ). As a consequence 
of the increased cooling efficiency (due to the increase of fjj 2 up to 1) the 
fragmentation process goes on until the Jeans mass reaches ~ 0.1 M Q or less, 
which is lower than the minimum mass for nuclear burning. 

As a result of the above picture, dark clusters with MHO of mass ~ 0.1M© 
or less would form in the outer part of the galactic halo. However, we don't 
expect the fragmentation process to be able to convert the whole gas mass 



contained in a PGC cloud into MHO (see, e.g., [P|)- Thus, we expect the 



remaining gas to form self-gravitating H 2 molecular clouds since in the 
absence of strong stellar winds (which in stellar globular clusters do eject the 
gas) the surviving gas remains gravitationally bound in the dark cluster. The 
possibility that the gas is diffuse in the dark cluster is excluded due to its 
high virial temperature (~ 10 4 — 10 5 K) that would make the gas observable 
at 21 cm. In addition, the gas cannot have diffused in the whole galactic 
halo because it would have been heated by the gravitational field to a virial 
temperature ~ 10 7 K and therefore would be observable in the X-ray band 
(for which stringent upper limits are available). The further possibility that 
the gas entirely collapsed into the disk is also excluded because then its mass 
would be of the order of the inferred dark halo mass. Finally, we expect the 
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ratio between the gas in the form of H2 molecular clouds and of MHO to 
be a function of the galactocentric distance which determines the ambient 
conditions such as incoming UV fluxes and collision rates among PGC clouds. 

A few comments are in order. The UV flux incoming in PGC clouds is 
crucial for determining fu 2 from which their subsequent evolution strongly 
depends. Since the formation of stellar globular clusters requires a sufficiently 
high UV flux, they can mainly form up to a certain galactocentric distance 
Rcrit which we estimate to be ~ 10 — 20 kpc |T2|. Beyond this distance the 
evolution of the PGC clouds gives rise to the formation of dark clusters of 
MHO and/or H2 molecular clouds. 

A further question which naturally arises is whether dark clusters are 
stable within the lifetime t g of the Galaxy. Standard calculations for the 
evaporation rate due to two-body encounters gives t evap ~ 100 t re i, where t re i 
is the usual two-body relaxation time. The condition t evap > t g requires a 



constituent mass m < 10 M [14J]. Another mechanism which could destroy 
dark clusters are collisions among themselves. We find that clusters are 
disrupted if they are located within a certain galactocentric distance Rdi S ~ 
10 kpc ||12|| . From these considerations we conclude that dark clusters of 
MHO and/or H 2 molecular clouds can still be present today at distances 
larger than Rdis- 

Let us briefly discuss the possible signatures of the above picture. Besides 
detection of microlensing events - which has in fact been our main motivation 



- we should mention that Maoz |22| has recently considered the possibility 
to infer from microlensing observations whether MHO are clustered or not. 
Already a relatively small number of events would be sufficient to exclude 
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this possibility, while to confirm it more events are needed. 

A further signature of the presence of H 2 molecular clouds in the galactic 
halo should be a 7-ray flux produced through interaction with high energy 
cosmic ray protons. Cosmic rays scatter on protons in H2 molecules produc- 
ing 7r 0, s which subsequently decay into 7's. 

As a matter of fact, an essential ingredient is the knowledge of the cosmic 
ray flux in the halo. Unfortunately, this quantity is unknown and the only 
information comes from theoretical estimates |23 . More precisely, from the 



mass-loss rate of a typical galaxy, one can infer the total cosmic ray flux 
in the halo, which turns out to be F ~ 10 41 A~ al erg cmT 2 s _1 , where 
Ag a i is the surface of the galactic disk. Taking R ga i ~ 10 kpc, we get F ~ 
1.1 x 10 -4 erg cm' 2 s' 1 . However, this information is not sufficient to carry 
out our calculations, since we need the energy distribution of the cosmic rays. 
We assume the same energy dependence as measured on Earth and then scale 
the overall density in such a way that the integrated energy flux agrees with 
the above value. Moreover, we assume that the cosmic ray density scales as 
1/R 2 for large R (R is the galactocentric distance). The measured primary 
cosmic ray flux on the Earth is 

^ m . 2 / E \ ~ 2 ' 7 particles . , 

and then the corresponding integrated energy flux is ~ 5.7 xlO -2 erg cmT 2 s~ x 
(integration range: 1 GeV < E < 10 6 GeV). Hence, & C r(E) ~ 1.9 x 
10 -3 §cri an d then we get 

* cr {E,K)~* C rW£±^ , (3) 
where a ~ 5 kpc is the halo core radius and Rqc ~ 8-5 kpc is our distance 



from the galactic center. The source function g 7 (r) giving the photon number 
density at distance r from the Earth is 

5 » = £ fdE r dE, ^cn(E r ,R{r)) ^AMl <K=fA (£ y , 

(4) 

where a^^{E^) is the cross section for the reaction pp — > nn (n is the 7r° 
multiplicity), n 7 (-E p ) is the photon multiplicity, -R(r) is the galactocentric 
distance as a function of r, while p^ 2 (i?(r)) is the fraction of dark matter in 
form of H2 molecular clouds (for which we assume the same dependence on 
R as in eq. (|]) ) . Actually, the cosmic ray protons in the halo which originate 
from the galactic disk are mainly directed outwards. This fact implies that 
also the induced photons will leave the Galaxy. However, the presence of 
magnetic fields in the halo might give rise to a temporary confinement of the 
cosmic ray protons similarly to what happens in the disk. In addition, there 
could also be sources of cosmic ray protons located in the halo itself, as for 
instance isolated or binary pulsars in globular clusters. Unfortunately, we are 
unable to give a quantitative estimate of the above effects, so that we take 
them into account by introducing an efficiency factor e, which could be rather 
small. In this way the 7-ray photon flux reaching the Earth is obtained by 
multiplying g 7 (r) by e/Airr 2 and integrating the resulting quantity over the 
cloud volume along the line of sight. The best chance, if any, to detect the 
7-rays in question is provided by observations at high galactic latitude, and 
so we take 9 = 90°. Accordingly, we find for the 7-ray flux |T2] 

$ 7 (90°) * e 1.7 x 10- 6 Ph f° ns . (5) 

cm 2 s sr 

The inferred upper bound for 7-rays in the 0.8-6 GeV range for high galactic 
latitude is 3 x 10~ 7 photons cm~ 2 s _1 sr' 1 Hence, we see from eq. (0) 



that the presence of H 2 molecular clouds does not led at present to any 
contradiction with the upper bound provided e < 10 _1 . 

In conclusion, we have outlined a scenario for a baryonic dark halo in 
which the formation of MHO and/or H 2 molecular clouds in the outermost 
part of the halo naturally arises, because in a quiet ambient the Jeans mass 
can reach values as low as 10~ 2 M . What is crucial in discriminating the 
evolution of PGC clouds towards stellar globular clusters or dark clusters 
are the decreasing collision rate and UV flux with increasing galactocentric 
distances. The most promising way to detect dark clusters of MHO is through 
correlation effects in microlensing observations, as we expect the dark clusters 
not to have been disrupted up to now. A much more difficult task is the 
detection of H 2 molecular clouds. A possible signature of such clouds would 
be their 7-ray flux induced by halo cosmic ray protons. Our calculation gives 
only an upper bound that is not in conflict with present detection limits. 
At any rate, the fact that our scenario naturally explains the formation of 
MHO at large galactocentric distances strongly supports - we believe - the 
idea that galactic dark matter is mostly baryonic, after all! 

We would like to thank B. Bertotti, F. Palla and D. Pfenniger for useful 
discussions. 
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